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Abstract: This paper is devoted to the adaptation of the method developed in Dolbeault et al. 2011 2009 to a 



Fokker-Planck equation for fiber lay-down which has been studied in Bonilla et al. 2007/08 



Gotz et al. 



2007 



Exponential convergence towards a unique stationary state is proved in a norm which is equivalent to a weighted L 
norm. The method is based on a micro / macro decomposition which is well adapted to the diffusion limit regime. 
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1 Introduction 

The understanding of the shapes generated by the lay-down of flexible fibers onto a conveyor belt is of great 
interest in the production process of nonwovens textiles that find their applications, e.g., in composite materials 



like filters, textile and hygiene industry. In Gotz et al. 2007 a stochastic model for the fiber lay-down process, 
i.e. for the generation of a fiber web on a conveyor belt has been presented. Taking into account the fiber motion 
under the influence of turbulence, the process can be described by a system of stochastic differential equations. 
We shall focus on a very simple situation, which does not take into account the movement of the belt, and gives 
rise to a Fokker-Planck equation. Some numerical results will also be considered at the end of this paper. 

An important criterion for the quality of the web and the resulting nonwoven material is how the solution 
converges to equilibrium. In particular, the speed of convergence to the stationary solution is important. The 
faster this convergence is, the more uniform the produced textile will be. From a technological point of view, 
process parameters should be adjusted such that the speed of convergence to equilibrium is optimal. 

The trend to equilibrium for solutions of kinetic equations has been investigated in many papers using 



entropy methods; see for example Desvillettes and Villani 2001 , Villani 2009 . A simplified approach has been 



suggested in Dolbeault et al. 2011 2009 . This trend to equilibrium for the Fokker-Planck equation for fiber 



lay-down under consideration in this paper has already been investigated using Dirichlet forms and operator 



semi-group techniques in Grothaus and Klar 2008 ; an ergodic theorem and explicit rates of convergence have 

an exp onential rate towards a uniq ue 

to 



at 



been established. In the present paper we prove the convergence 

stationary state in a weighted L 2 norm by adapting the method developed in Dolbeault et al. 
the setting of non-moving belts. 



2011 



2009 



2 The model and main results 



In the melt-spinning process of nonwoven textiles, hundreds of individual endless fibers obtained by the 
continuous extrusion through nozzles of a melted polymer are stretched and entangled by highly turbulent 



air flows to finally form a web on a conveyor belt (see Gotz et al. 2007 for more details). We describe the 



motion of an individual fiber, neglecting interactions with the others. 
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An arclength parametrization of the laid down fiber in a coordinate system following the conveyor belt 
is given by x (t) e M 2 , t > 0. The tangent vector is denoted by dxo(t)/dt = r{a{t)) with r(a) — (cosa,sina), 
a€ S 1 = M/27rZ. Since the lay-down process is assumed to happen at the constant normalized speed 1 (equal 
to the spinning speed), xo(t) can also be interpreted as the position of the lay-down point at time t. If the 
conveyor belt moves with velocity ne\, the history of the lay-down point in the laboratory frame (as opposed 
to the conveyor belt frame) is given by x(t) — xo(t) + tKe\, i.e., 

^ =r(a)+Kex. (1) 

It is a natural restriction that the speed of the conveyor belt cannot exceed the lay-down speed: < k < 1, since 
otherwise a stationary lay-down point would be impossible. The lay-down process can now be determined by 
prescribing the dynamics of the angle a(t), decribed as a stochastic process. It is driven by a deterministic force 
trying to move the lay-down point towards the equilibrium position x — and by a Brownian motion modeling 
the effect of the turbulent air flow: 

da = -T ± (a)-W(x)dt + AdW, (2) 

where W denotes a one-dimensional Wiener process, A > measures its strength relative to the deterministic 
forcing, r = dr /da — (— sin a, cos a), and V(x) is a potential such that e~ v is integrable with the normalisation 
/ R2 e~ v dx = 1 and W(Q) = 0. 

The system Q-Q defines a stochastic process on M 2 x S 1 . The corresponding probability density f(t, x, a) 
satisfies the Fokker-Planck equation 

d t f + (r + K ex) ■ V x f - d^r 1 - -V x Vf + Dd a f)=0, 

(with the diffusivity D = A 2 /2) which will be the object of our study. The analysis of the long time behaviour 
is considerably simplified in the case of a nonmoving conveyor belt: 

k = 0. 

We shall assume that this assumption holds true from now on. In this case the Fokker-Planck equation is written 
as an abstract ODE 

d t f + Tf = DLf, (3) 

with T/ = r • V x / — d^r 1 - ■ V x Vf) and If = d 2 J. 

It is easily seen that F(x, a) = e~ v ( x ) is an equilibrium solution of lying in the intersection of the null 
spaces of T and L: TF = LF = 0. 

A convenient functional analytic setting is introduced by the scalar product 

(f,9)-=f fgdu, dfj,(x,a) := ^f" °[ , dv(a) := ^ , 
y R2xS1 F{x,a) 2tt 

and by the associated norm ||/|| 2 = (/,/)• On the space L 2 (M 2 x S 1 ,dfi), the operator T is skew symmetric, 
and the operator L is symmetric and negative semi-definite. Thus, we have 

j t ^p^=D{Lf lf ) = -D\\d a .ff. (4) 

This identity reveals the main difficulty in proving convergence to equilibrium. The decay to equilibrium seems 
to stop, as soon as / is in the null space of L consisting of all a-independent distributions. On the other hand, 
the decay equation Q does not make use of the action of the operator T and, in particular, of the fact that the 
equilibrium F is the unique probability density in Af(T) D A/"(L). Any a- independent distribution function / is 
indeed unstable under the action of T, unless f = F. Hence, convergence to the equilibrium can be expected 



and will be proven to be exponential. This is a so-called hypocoercivity result as defined in Villani 2009 1. A 
recently developed approach Dolbeault et al. 2011 2009 for proving hypocoercivity in the abstract setting ^ 



will be applied with a special emphasis on the behaviour of the decay rate as D — > and D — > oo. It requires 



assumptions on the potential V, which have already been used in Dolbeault et al. 2011 



(HI) Regularity: V £ W^{R 2 ). 
(H2) Normalization: J K2 e~ v dx = 1. 
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(H3) Spectral gap condition: there exists a positive constant A such that 



/ 



\VxU\- 



dx > A 



2 -V 

u e 



dx 



for any u G H 1 (e v dx) such that Jj 



u e 



dx = 0. 

(H4) Pointwise condition: there exists c\ > such that 
\V 2 x V{x)\ < ci (1 + \V x V{x)\) for any a; € R 2 . 

Roughly speaking, (H2) and (H3) require a sufficiently strong growth of V(x) as \x\ — > oo, whereas (H4) puts a 
limitation on the growth behavior. This leaves room, however, for a large class of confining potentials including 
V(x) = (l + \x\ 2 )P, 0> 1/2. 



In 



Dolbeault et al 



2011 



the additional pointwise condition A x V(x) < | IV^I^iz;)! + cq with cq > and 
€ (0, 1) has been required. This is, however, a consequence of (H4) by A x V(x) < ^2 |V 2 V»| and the Young 
inequality V2ci \V X V\ < \V x V\ 2 /<i + 2 c 2 , with = 1/2 and c = c x + 2 c 2 . 

Theorem 1. Let f G L 2 (R 2 x S 1 , dfi) and let (HI )-(H4) hold. Then, for every n > 0, the solution of ^ subject 
to the initial condition f(t = 0) = f satisfies 



\\f(t)-F\\ < (1 + V )\\f -F\\e 



-At 



with A = 



n 



CiD 



1 + n 1 + C 2 D 2 

where C\ and C2 are two positive constants which depend only on the potential V . 
As a consequence, we have 

A = 0(D) asD-s-0 and A = 0(D^ 1 ) as D -> 00 . 



Both results are sharp, as can be seen from the toy problem in Dolbeault et al. 2011 . As D — > 0, dissipation 
is provided by the 0(D) right hand side of (j3|), which dominates A. On the other hand, the dynamics for large 
D can be described by a macroscopic limit: see |Bonilla et al. |2007/08| . When D —> 00, the correct time scale is 
t = O(D) and corresponds to a parabolic scaling, and therefore A = 0(D~ 1 ) had to be expected. 

However our method is not sharp in the sense that we cannot expect to obtain the optimal coefficients in 
the limiting cases above. Again this can be seen from the toy problem in Dolbeault et al. 2011 , where the 
spectral gap can be explicitly computed. 



3 Proof of Theorem Q] 

3.1 The modified entropy 



We introduce the deviation g := f — F, satisfying ^ subject to g(t = 0) = /o — F. Following Dolbeault et al. 



2011 , we denote the orthogonal projection to A/"(L) by 



n<? ■= Pg= / gdv. 
Js 1 

Note that J ffi2 p g dx = holds. In the following, we shall also need 

Tn 9 = r ■ (V xPg + p g V x V) = r ■ e- v V x (e v p g ) 

with the consequence 

nm = , 



(5) 



which is essential for the applicability of the method of Dolbeault et al. [2011 
limit (corresponding to D — > 00) in Q is diffusive: see Bonilla et al. [2007/08 
With the help of the operator 

a = (1 + (Tn)*Tn) _1 (Tn)* 

and an appropriately chosen e > 0, the modified entropy functional is defined by 



(G) 

It implies that the macroscopic 



H[ff] :=\\\g\\ 2 + e(Ag,g) 
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By Dolbeault et al. 2011 Lemma 1], (|6| implies 



|A 5 ||<i||(i-n) 5 || 



Hence the modified entropy functional is bounded from above and below by the square of the norm for any 
£ G (0, 1). More precisely we have: 

7 Sll 2 < H[ff] < 



1 1 M ^ u „ ^ „ ^ 



2 

A straightforward computation gives 



±H[g] = -D[g], (7) 



with the entropy dissipation functional 

D[g] = -D (Lg, g) + e (ATIlg, g) + e (AT(1 - II) 5 , .9) - e (TAg, <?) - e D (ALg, g) . (8) 

3.2 Microscopic and macroscopic coercivity 

The first term on the right hand side of ^ has already been computed in Q. With dv = da/ (2tt), the Poincare 
inequality on S , 



's 1 Js 
implies the microscopic coercivity property 



/ \d a g\ 2 dv > I (g - J ' g dv) 2 dv 
Js 1 Js 1 



<L 5 , 5 > > ||(i-n) ff r . (9) 



The operator ATI! = (1 + (Til)* Til)- 1 (Til)* Til shares its spectral decomposition with (Tn)*TII. For the latter 
we have, using ^ 

((Tn)*Tn 5 , g ) = ||Tn.g|| 2 = i / e - y \V x u g \ 2 dxdv, 
with u g = e v p g . The spectral gap condition (H3) implies the macroscopic coercivity property 

((Tn)*Tn.g,. 9 )>|||p s || 2 , 

leading to 

(ATn g ,.g)>^_||n.9|| 2 . (io) 

By (9) and the sum of the first two terms in the entropy dissipation (fST) is coercive. This will also be 

sufficient for controlling the remaining three terms, if the operators AT, TA, and AL are bounded, for e > 0, 
small enough. 

3.3 Boundedness of auxiliary operators 

By |Dolbeault et al. 2011 Lemma 1], we know that 



||TA 5 |j<||(i-n) 5 ||. (ii) 

The computation 

(TUyig = -UTLg = -V x ■ U(rd 2 a g) = V x ■ U(t g) = -(TU)*g 
shows that AL = —A and, thus, 

||AL fl ||<i||(l-n) 5 ||. (12) 
The most elaborate part of the analysis is to prove the boundedness of AT. Following the approach of Dolbeault| 



et al. 2011 , we consider its adjoint 

(at)* = -T 2 n(i + (Tn)*Tn) _1 
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For a given g £ L 2 (R 2 x S 1 , dp), we introduce h = (l + (TIT)* Til) 1 g which, after solving for g and applying IT, 



becomes 



Pg = Ph - ^ 2 Ph = e v u h - - V x ■ (e v V x u h ) (13) 



with Uh — e v ph- A straightforward computation gives 

(AT)*g = -T 2 Ph = e- v [(t ■ V,) 2 u h - (r x • V X V) (r x • V*^)] 

and, as a consequence, 

||(AT)*. 9 || < IIVfVll^^^v^ + lllV^llV^HI^^^v^ . 

Therefore, in order to prove the boundedness of (AT)* (and, thus, of AT), we need to prove the boundedness of the 
right hand side in terms of \\p g \\ for the solution Uh of the elliptic equation ( 13 ). This L 2 — > H 2 (with weight e~ v ) 
elliptic regularity result has been derived in Dolbeault et al. [201 1| under the assumptions (H1)-(H4) (see 
Proposition 5 for the first term and Lemma 8 for the second). Collecting these results gives ||(AT)*g|| < Cy \g\ 
and therefore 

||AT(i-n). 9 || <cv||(i-n) fl ||, (14) 

where Cy depends only on the potential. 
3.4 Hypocoercivity 

Inserting estimates (|9])-( 12 ) and ( 14 1 in Q gives 

D[g] > D\\(l -n). 9 || 2 + ^||n 5 || 2 -e(Cv + l+f) ||(1 - n) 5 || || 5 || 

> (d - s (c v + 1 + f )) ||(i - n). 9 || 2 + ^_ ||n 9 || 2 - e (cv + i + f) 11(1- m\ \m (15) 

> {D - e {Cy + 1 + f ) (1 + i)) ||(1 - IT). 9 || 2 +£ (^ K -I (C V + 1 + f )) ||II 5 || 2 , 

for an arbitrary S > 0. This shows already that coercivity can be achieved by first choosing S and then e, both 
small enough. With the choice 

6= A 



(2 + A) (Cy + 1 + D/2) ' 



the coefficients on the right hand side of (151 can be written as D — e r(D) and e s with 
r(D) := (2A + (2 + A) (Cy + 1 + f )) (Cy + 1 + f ) , 



2A v v y v " 2 " y v 21 ' 2(2 + A) ' 

Then the optimal choice of e, considering the form of the coefficients, would be 

However, we also have to guarantee e < 1 for the definiteness of H[g) and actually, even stronger, y±| < (1 + r)) 2 
will be needed below, which can be guaranteed by the requirement e < Moreover the two conditions are 
equivalent at first order for r/ > 0, small. These considerations lead to the choice 

e = , with e m ax = max { 1 , max e(D) ) , 

1 + Tj Zmax D>0 > 

which is finite because of e(0) = e(oo) = 0. With this choice, 

with appropriately chosen constants C\, Cy. > 0, depending only on A and Cy and, thus, only on the potential V. 
The estimate 

D[ 5 ]>2A|| 3 || 2 > I ^H[. 9 ]>2AH[. 9 ] 
follows. Using this in Q and the Gronwall lemma imply 

H [/(<)-*! < H[ib-F]e- 2At . 

Finally we obtain for the norm 

\\f(t)-F\\ 2 < Y ^- £ H{f(t)~F}< I ^- £ H{f Q -F]e- 2 ^ 

< ^ ||/o - ^ll 2 e- 2At < (1 + vf ll/o - ^ll 2 e- 2At , 

1 — £ 

which completes the proof of Theorem [T] 
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Figure 1: Plot of t^> log \\f(t, ■, •) - F|| 2 for A = 0.5, 3, 6. 



4 Concluding remarks 

4.1 Numerical investigations 

It is interesting to compare the rates predicted by the above results, which are only upper bounds, with numerical 
rates of convergence. We use a classical Monte-Carlo method with an Euler-Maruyama discretization scheme 
for all computations. A numerical investigation of the equations using a semi-Lagrangian method can be found 
in |Klar et al] |2009 . 

The exponential decay of the L 2 -difference to the stationary solution is observed in Figure [TJ In Figure [2] the 
decay rates A have been obtained from the above simulations for various values of A using a least square fit. The 
rate given by TheoremJTJ i.e. A ~ TTcfl? 7 ' ^ = A 2 /2, fits qualitatively very well the curve obtained in Figure^ 
when A is away from CTTIn particular, values of A with an optimal rate of convergence can be determined from 
the numerical as well as the analytical results. 



4.2 Perspectives 

1. Models where stationary solutions are not known explicitly. The model considered in this paper can be 
extended in different directions, for instance by taking into account the movement of the belt, or by 
models where fibers have smoother trajectories than the ones considered in this paper, see |Herty et al.| 
2009] , [Bomlla et al.| [2007/08] . In these cases the stationary solutions are not always known explicitly. The 



application of the entropy method presented above is then an open problem. 
3-D models. To model fluid flow through a fiber web, the model has to be extended to three dimensional 
situations, see Klar et al. 2012 . For such a model exponential convergence to equilibrium, at least for the 



case k = 0, can be proven with the same methods as in this paper. 
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